EQUIDISTANCE OF THE COMPLEX 2-DIM ANHARMONIC 
OSCILLATOR SPECTRUM: EXACT SOLUTION 
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We study a class of quantum two-dimensional models with complex poten- 
' r~ | '■ 

tials of specific form. They can be considered as the generalization of a 
| recently studied model with quadratic interaction not amenable to conven- 

tional separation of variables. In the present case, the property of shape 
invariance provides the equidistant form of the spectrum and the algorithm 
to construct eigenfunctions analytically. It is shown that the Hamiltonian is 
non-diagonalizable, and the resolution of identity must include also the cor- 
responding associated functions. In the specific case of anharmonic second- 



o 

[ plus-fourth order interaction, expressions for the wave functions and as- 

o 

CN ' sociated functions are constructed explicitly for the lowest levels, and the 



recursive algorithm to produce higher level wave functions is given. 
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Introduction. 



Considerable attention has been paid recently to the quantum models with complex poten- 
tials, and in particular, to the models with real spectra pp. Although analysis of such systems 
is fraught with many difficulties [2], [3], it seems to be very useful for modern development 
both of generalized Quantum Mechanics and Quantum Field Theory [I]. In particular, a 
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modified scalar product which provides unitary evolution must be built, and diagonalization 
of the Hamiltonians has to be controlled [5j. Very interesting results were obtained for a 
wide class of models with unbroken PT— invariance [1], [6] - [8]. The important notion of 
the pseudo-Hermiticity was introduced [9] - [TTj : 



with r] a Hermitian invertible operator. It allowed to define a class of non-Hermitian systems 
with physically acceptable properties of energy spectra. Systematically pseudo-Hermiticity 
has been presented by A.Mostafazadeh [9] (see also [10] - [12] ). The Hilbert space for such 
kind of systems is described in terms of a biorthogonal basis, which consists of the eigenstates 
|\I/ n ) and \^ n ) of H and W, correspondingly. The essential peculiarity of non-Hermitian 
models is known: some of them are naturally described by non-diagonalizable Hamiltonians 
[5]. For such systems, the usual biorthogonal basis does not provide a complete basis in 
Hilbert space, and one has to add the so-called associated functions to complete the basis. 
Then, the Hamiltonian is block- diagonal with some number of Jordan blocks of standard 
structure on its diagonal. 

It must be noticed, that during last years almost all papers in this line were limited 
to one-dimensional Quantum Mechanics, where the problems have more chance to be fully 
solved. Nevertheless, some investigations of two-dimensional systems were performed as 
well [13] - [15] . Recently, a new non-trivial two-dimensional model with complex potential 
was studied in detail in paper [16]. The potential of this model had the simple form of 
second order polynomial in X\,X2- Usually models of such kind are solved by means of linear 
transformation of coordinates with subsequent separation of variables (see for example, [H]), 
but for two values of coupling constants, the Hamiltonians are not amenable to separation 
of variables. Just such model was studied in [16]. Exact solution was possible essentially 
due to shape invariance property [17] - [23] originated from the Supersymmetrical Quantum 
Mechanics approach [21], |25j . For the model [16], the energy spectrum was proven to 
be equidistant, the same as for usual real harmonic oscillator. The corresponding wave 
functions were found analytically as well. From the properties of wave functions it became 
clear that the Hamiltonian is non-diagonalizable, and its structure was analyzed in detail. 
The corresponding associated functions were also built analytically and their properties were 
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studied. 

In the present paper, the generalization of the model [16] is given: new exactly solv- 
able models are built with the full spectrum and the wave functions analytically calculated. 
In Section 2, we consider the wide class of complex two-dimensional potentials which in- 
clude two terms: an arbitrary polynomial function of z = X\ — ix 2 in addition to real 
two-dimensional isotropic harmonic oscillator. Somehow unexpectedly many properties of 
this system will turn out to be analogous to those of complex quadratic potential of [IB] . 
Namely, the equidistance of spectrum and analytical expressions for wave functions follow 
directly from the shape invariance. In this sense, the situation is very different from that of 
the Hermitian case, where the spectrum of anharmonic oscillator (which does not obey any 
shape invariance property) can be very far from the equidistant. In Section 3, the differ- 
ent properties of the system are studied, such as its pseudo-Hermiticity and modified scalar 
product. The self-orthogonality of all excited wave functions signals non-diagonalizability of 
the Hamiltonian with all related features. The analysis leads to the specific algorithm for 
construction of associated functions for all Jordan blocks for the particular case of specific 
quartic anharmonic form of interaction. In Section 4, two lowest Jordan cells are studied 
explicitly, all associated functions are constructed, and scalar products are shown to satisfy 
the standard rules for non-diagonalizable pseudo-Hermitian Hamiltonians. 

2. Description of the model. 

Let us start from the two-dimensional Hamiltonian with complex interaction of the form: 

H = -A (2) + V(x) = -Ad z d- Z + X 2 zz + 2XzF'(z); (2) 
z = x\ + ix 2 ; z = x\ — ix 2 ; F'{z) = d z F(z) 

where F(z) is supposed to be a polynomial. This potential is the direct generalization of 
interaction in [16] . where F' was proportional to z with a coefficient such that the con- 
ventional separation of variables is impossible. Also the present case (T5]) is not amenable to 
separation of variables [26], [27], and the spectrum of the model again can be found by means 
of shape invariance [17], [19] - [21]. This property was introduced [17] in the framework of 
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Supersymmetrical Quantum Mechanics as 

H {2) (x; a) = H {1) (x; a) - K(a) (3) 

between two partner Hamiltonians which are intertwined by supercharge operators : 

H (1 \x; a)Q + (x) = Q + (x)H^{x; a); Q-(x)H (1 \x; a) = H^\x] a)Q~(x). (4) 

Above, a and a = a{a) are some parameters which define coupling constants of H^ l \ H^ 2 \ 
and lZ(a) does not depend on x. Shape invariance allows to find algebraically the spectrum 
and wave functions of Hamiltonians - completely in one- dimensional case p2], [18], and 
partially - in two-dimensional situation [19], [20], [22J. 

The simplest realization of one- dimensional shape invariance with a = a, i.e. with 
H (2) (x,a) = H {1) (x;a) + 2A, was considered in [21] . For that case, both systems were 
shown to have oscillator-like (equidistant) spectra. The same kind of shape invariance and 
of spectra appeared in two-dimensional systems with complex quadratic interaction [16] . 
The basic origin for this property of spectra are the following - oscillator-like - commutation 
relations for H : 

HA + = A + (H + 2A); HA" = A~{H — 2A). (5) 

These relations are the particular case of (JIJ), where A play the role of supercharges Q ± , 
2A - the role of IZ(a), and H, H + 2 A - the role of partner Hamiltonians H^'^ 2 \ In general, 
an arbitrary system with Hamiltonian obeying commutation relations fl3]) has oscillator-like 
equidistant spectrum [2T] , maybe on a finite interval of energy: its lowest level is defined 
by the zero mode of A~ and the highest level - by the zero mode of the operator A + . Since 
the intertwined Hamiltonians in (J3]) coincide up to a shift, they are called self-isospectral. 
Thus, both for general one-dimensional [21] and for specific two-dimensional [16] cases, self- 
isospectrality leads to an equidistant (oscillator-like) character of the spectrum. 

Let us consider now in detail the system with Hamiltonian (j2J) . Due to absence of singu- 
larities and similarly to the model [16], no "fall to the center" phenomena [28] is possible here 
(see the formal proof in Appendix of [16]). The normalizable bound state wave functions will 
be exponentially decreasing at infinity, having no singularities. The corresponding spectrum 
is bounded from below, and the ground state with energy E will be denoted as \E , ,o( : ?) (h 
will be clear below, why we use two indices for enumeration of \1/). 
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First of all, the Hamiltonian (TSl) satisfies the commutation relations (El) with first order 



operators A ± of the form: 



A ± = d z T \z. (6) 



Similarly to [21], [16], it is necessary to define the zero modes of the intertwining operator 
A- : 

A-* n , (z,z) = Q. (7) 
Looking for the zero modes in the form: 

V n> a(z, z) = exp (~zz - F(zU * n , , (8) 

one easily obtains that ^ n ,o depend only on z. The zero mode ty nt o(z, z) to be simultaneously 
eigenfunction of H with energy E n has to satisfy the simple equation for ty nj o(z): 



2A 

Its solutions are: 



zd- z + 1 



V n ,o(z) = E n y n> o{z). 



E n —2X 



%n,o( Z ) = °nfiZ 2X , (9) 



where c ni o are constants. Since the wave functions must be the single-valued functions on a 
plane, i.e. to be 2ir— periodic in polar angle ip, the allowed energy spectrum of the model is: 

E n = 2\(n+1); n = 0,1,2,.... (10) 

The ground state corresponds to n = 0, and excited states - to n = 1, 2, .... One can argue 
that no additional bound states are possible besides (flO]) . Indeed, supposing that such 
additional level with energy E exists, one may act several times on its wave function by the 
operators A~ . According to the second relation in (jSJ), one will obtain wave functions with 
energies E — 2Am; m = 1, 2, .... This series can be cut only if the zero mode of A~ appears 
on some stage, but all zero modes are described by (jHJ), © above. Therefore, the ground 
state of the model is: 

%,o(z,z) = co t oexp(--zz-F(z)J; E = 2A. (11) 
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One can check by direct calculation that wave functions of excited states given by (jHJ), © 
can be also written as: 

V nfi (z, z) = c nfi z n exp (~zz - F(zfj = (A+) n ^ 0>0 (z, z), E n = 2\{n + 1); n = 1, 2, 

(12) 

with c„ i0 = (-A/2) n c , . 

Thus, the whole spectrum ( FlUj) and the corresponding wave functions ty n>0 (z,z) of the 
system (j2j) are known analytically. It should be noted that the spectrum (ITU]) depends neither 
on the structure of function nor on its coupling constants (coefficients b and u below 
Eqs. (|27p and fl3Tj) ). This is a general consequence of relations (jSJ). A further investigation of 
properties of the wave functions will be performed in the next Section. 

3. Non-diagonalizability. 

It is known [2], [5] (see also [TB]) that the self-consistent formulation of Quantum Mechan- 
ics with pseudo-Hermitian (TjQ) Hamiltonians includes a suitable modification of the scalar 
product and resolution of identity. Namely, a new scalar product has to be used: 

= J (i/T¥)S, (13) 

so that the pseudo-Hermitian H becomes Hermitian when equipped with fll3p . In our present 
case (ED, one may define rj = P 2 , the reflection x 2 -H- — x 2 , and therefore, the scalar product 
(T3~3|) is simply an integral over the product J instead of the J in the ordinary 
Quantum Mechanics, where rj = 1. 

First of all, we have to check the properties of all wave functions ( fl2l) in the framework of 
a new scalar product. Since below we will need the results of integration of some functions 
over z, z, the list of such integrals is given in Appendix. In particular, the norm of states 
^ n ,o vanishes for n > 1 : 

r 2 r 

\2j2^ n,0 / -2n 



((^ n ,o|^n,o)) = / {^nYd 2 x = / z 2n exp [-Xzz - 2F(z)]dzdz 



exp(-2F(0))5 n0 (14) 



A 

due to expansion of exp (2F(z)) in powers of z (see the list of integrals in Appendix). 

6 



Thus, the usual resolution of identity in terms of the so-called biorthogonal basis [3], [9], 
[29], [15], [16] l^n)) l^n) does not hold, and the Hamiltonian is non-diagonalizable. Details 
can be found in [5] in a one-dimensional case, and in [16] in a two-dimensional context: it 
is necessary to build the so-called associated functions which participate to the resolution of 
identity and complete the basis. 

For such systems, each self-orthogonal wave function \I/ ni o, n > 1 with zero norm must 
be accompanied with a set of p n — 1 associated functions \l/ ni fc, k = 1, 2, ...,p n — 1, where p n 
is called the dimension of Jordan cell. This situation (c.f. Eq. (fl5|) below) is very different 
from the familiar one with degeneracy of energy level of Hermitian Hamiltonian. It should 
be clear now, why notations with two indices of wave functions were introduced above. By 
definition, these functions obey: 

(H-E n )V ntk = V n , k -i; k = l,2,..., Pn -l, (15) 

where all functions are supposed to be normalizable. Each self-orthogonal eigenfunction 
^n,o, n = 1,2,... is supposed to be accompanied by only one set of associated functions 
*„ >fc , k = 1,2, ...,p n - 1. 

Similarly to the scheme of the previous paragraph, the partner eigenfunctions \l/ ni o also are 
accompanied by their associated functions ^ n ,fc, k — 1,2, ...,p n — l. The following numeration 
for the functions \1/ is convenient: 

*n^-*-i = Kk k = 0,1,2,. ..p n -l. (16) 

With these notations, according to the general formalism which was illustrated in detail for 
some one-dimensional models [5], the scalar product in the extended biorthogonal basis is: 

«»„!*-,» = <*■,!*-,> = /M!)M!)A = <") 

k = 0, 1, ...,p n - 1; / = 0, 1, ...,p m - 1. 
Correspondingly, the generalized decompositions become: 

OO Pn — 1 

n=0 fc=0 

CXD Pn — 1 OO Pn— 2 

^ = EE^|^))((^,Pn-fe-l|+EEl^,fc))((^,p„- fc - 2 |. (19) 

n=0 fc=0 n=0 fe=0 
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The Hamiltonian H is clearly non-diagonal, but block-diagonal. Each block - Jordan cell of 
standard form (see ([19]) ) - has dimensionality p n , and we assume the absence of degeneracy: 
each eigenvalue corresponds to one Jordan cell. 

The crucial property of eigenf unctions ^ n) o - their self-orthogonality for n > 1 - was 
already displayded in (TH|) . To complete the construction of the Jordan cell, it is necessary 
to find the corresponding associated functions with properties listed above. 

The scalar products ( IT71) vanish for different energy levels E n , E m , i.e. ({9 n ,k\^rn,i)) — 
for n m. For the proof, the pseudo-Hermiticity of H is important. For example, similarly 
to the ordinary Quantum Mechanics, wave functions with different energies are orthogonal: 

= ((V ni0 \H\y mj0 )} - ((## n) o|# OT , }) = (E m - £ n )((# n)0 |# OT ,o}). 

Further, because of (fT5l) : 

= ((# n , |(F - E m )\$ m>1 }) - (((H - E m )y n>0 \ty mj i)) = 

= «#»,o|*m,o)) - (K ~ E m ) ((^ol^l)) , 

and the scalar products between wave functions and first associated functions for different 
E n , E m also vanish: 

((*n,o|*m,l)) = «*n,l|*tn,o)) = 0. 

The procedure can be continued leading to orthogonality of all functions with different n, m. 

Such indirect method does not allow to study the scalar products of associated functions 
with the same value n : it is necessary to build them explicitly. This task is more difficult 
than in the case of harmonic oscillator [16] . Let us transform the defining equations to the 
form which can provide some simple algorithm to construct associated functions. It follows 
directly from the commutation relations ([5]), that 

(A-) k (H - E n ) = (H - E n . k )(A-) k . 
Together with relation (fl5l) . for k — 1 it gives: 

(H-E n - 1 )A-V ntl = 0, 
and therefore, the wave function and its first associated are related: 

A~^ n ,i = a n)1 * n _i j0 ; a n>1 = Const. (20) 
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Iterating this procedure, one obtains: 

{A~) k m n)k = a njk ^ n _ kfi ; a nM = Const, (21) 
with the last relation of this sort: 

(A-) n * n , n = a n , n * ,o- (22) 
Functions in r.h.s. of (1201) - (|22|) are known from (Tl2|) . and the ansatz similar to ([8]) 

^„ ifc = exp (-^zz - F(z))^ njk 
is useful to look for all associated functions. Indeed, one can check that 
(A-) m exp {-\zz - F(z)) = exp {-\zz - F{z))d™, 
providing the following inhomogeneous differential equation: 

Its general solution is the polynomial in z : 

fe-i 

^n,fc - jTj Z Z + Z^ g nA Z ) Z ' y Z6 > 

but the functions gn\{z) are yet arbitrary functions of z. These functions have to be found 
by substitution of (123|) into relations ( IT5]) connecting fc with ty n> k-i- Due to: 

(H-E n )exp(~zz-F(z)) = 2exp(~zz-F(z)) \^-2d z d z +Xzd z + Xzd z + 2F'(z)d z -Xn 
the equations for g^ k {z) take the form: 

2^-2d z d- z + Xzd- Z + Xzd z + 2F\z)d z - Xn^j ( ^h^M z ^n-k + ^ = 

fc-2 

dn,k-\C n -k+lfi k-l-n-k+l 



-Z Z 



(k-l)\ 
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Considering fl2l]) as polynomial in z, one has trivial identity in highest power z k , but for the 
next powers, fl2^|) gives the system of coupled first order differential equations for unknown 
function g^\(z), which are written separately for % = k — 1 as: 



2A 



d n ,k-lC n -k+l,0 i , 4a n k c n-h ,0 

-Z 4 



{k-l)\ 



{k-l)\ 



(n - k)z~ 3 - F'{z)z 



(25) 



and for i — 0, 1, {k — 2) as: 



2A 



z l - n 9 { X{z) 



=*— n— 1 



(26) 



It is convenient to extract explicitly the quadratic term from the function F(z) : 

b 



F(z)="-z 2 + f(z). 



(27) 



Then, equation ( 125]) takes the form: 
2A 



a n,k-lCn~k+l,0 ~ 46a n k c n-k ,0 1 . 4a n fcC^-ft n 

—z 4 



(*-l)! 



(*-l)! 



(n-Jfe)z- 3 -/'(z)z- 2 
(28) 

Since the wave functions, and therefore the functions g^k~ l \ must be single- valued on a 
plane, we have to exclude the logarithmic term from the solution of fl28|) . This is possible 
by choosing the coefficient in the first term in r.h.s. vanishing: 



and the solution of (1251) is: 



X(k - 1)! 



a n,fc-lC n -fc+l,0 — ^ a n,k c n-kfii 



a (fe-l)_„_ fe+1 + (n _ + 2 ^n- fc+ l / f'^)z~ 2 dz 



(29) 



(30) 



where aj^^ are integration constants. 

Now, (!30|) and (|26|) can be used to find iteratively functions g^ k {z) for 2 = k — 2, i = k — 3, 
and so on. Thereby, the required associated functions ^ n> k will be obtained. Below, for 
simplicity we restrict ourselves to the typical quartic anharmonic interaction: 



m = \z\ 



(31) 
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which will clarify the subsequent construction without too cumbersome formulas. With this 
interaction, expressions ( 130]) give the first associated function which depends on the only 
function g^\{z) : 



&n,l c n-l,0 



a nt ic n -i fi zz n 1 + gf\{z) 



1 ra— 1 ^n— 2 n,i -n 2uzn+2 



exp y-\{Xzz + 6z 2 + W2 4 )j = (32) 



Analogously, the second associated function depends on two functions g : 



n.2 



r 2-n-2 , „(0) 



a n ,2C n -2,oz z + g"i{z) + g y ' 2 [z)z 



v:\\) i —-(Xzz + bz 2 + loz 4 



Here, g^\ is given by ( l30l) . but g£\ has to be calculated from ( |26l) with i = 0, k 
taking into account relation ( |29i) : 



(33) 
2 and 



,-nJO) 
z i/n,2 



2a n .20 



n,2C„_2,0 



A 2 



(n - l)o^\ + b)z' 3 + (n-2)(n- 3)z~ 5 + 



Again, to avoid the logarithmic term in g^ 2 the relation 

6u = b(afy - c£\) 

must be fulfilled. 

Thus ( |33l) provides the explicit expression for the second associated function: 

1 ,v2=n-2 



„(1) ^ _ X,.,2=3 



(34) 



(35) 



+ A^ 



*n,2(«) = a n , 2 c n _ 2 ,oS |« ^ n + 



(n - l)c£\ + b]z n - 2 + |(n - 2)(n - 3)^~ 4 + 2waJ l 1 ^ n+2 + 2w 2 ^ n+4 + 2/3 n)2 z* 



(36) 



^ | exp ^-§(Azz + bz 2 + ojz 4 )^ 



Higher associated functions can be calculated straightforwardly with an analogous procedure. 
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4. The lowest Jordan cells. 



For n = 1, 2, results ( 1321) and ( 1361) obtained in the previous Section are already sufficient to 
know the corresponding Jordan cells completely. For the first Jordan cell with dimensionality 
p n = n + 1 = 2, the wave function and associated function are: 

*i,o = ci j0 ^exp (—-(Xzz + bz 2 + uz 4 ) J (37) 



A ' 



exp ( — -(Azz + &z 2 + wz 4 ) ): .V|.j = Const. 



To be sure, it is necessary to check the scalar products inside the Jordan cell. We know 
already that the wave function is self-orthogonal. The scalar product ((^i,q|\Pi,i)) is: 



j # 1>0 *i,id 2 x = ^J y lfi y 1A dzdz = c 1}0 N hl J exp 



(Xzz + bz + ujz ) zz dzdz 



= ci, iVi,i exp (^-^-^j ^-J^J(*A c )|c=o = ^ci t0 N ltl , (38) 

where we used the fact that derivatives over b and A commute, and that all integrals In=o,m 
vanish due to relations in Appendix. Thus, choosing suitable normalization constants: 

«*i,o|*i,i» = 1, (39) 

and, what is important, it does not need fixing the constant a^ ]. 

For the last scalar product in the first Jordan cell, we have to calculate: 

= ^1,1 exp (-u-jjp) J ex P (-{>*z + fe2 )) (z 2 ~ j^f}zz) dzdz = 

= 5^«p(-a (-| + ? #^)^ i -^ = -^ + ^ (40) 

where again we used the integrals of Appendix. Here, the integration constant afl, which 
was not yet fixed, can be chosen as af'l = —6, so that the scalar product (1401) vanishes. 
Therefore, the first Jordan cell has dimensionality pi = 2, it includes the wave function and 
associated function: 

*i,o = ci, ^exp + bz 2 + uz 4 ^J (41) 

= jV^i (^Xz + bz- 2uz 3 J exp (^-^(Xzz + bz 2 + uz 4 ] 

12 



and all scalar products just correspond to (TTTj) . 

The situation with the next Jordan cell n = 2 with dimensionality p 2 = 3 can be studied 
analogously. 



^2,0 = c 2 ,o<2 2 exp ( --(Xzz + bz 2 + uz*] 



*2,1 = N 2A 



zz - -(1 + a ( 2 }z 2 + 2uz 4 ) 



exp ( —-(Xzz + bz 2 + uz 4 ] 



*2.2 = A 2.2 <( 



1 



+ 6) + 2wa 2 i 2 ^ 4 + 2wV + 2/3 2 , 2 z 2 



.,(1)^4 



-j(a^z + 2uz 3 



exp ( — -(Azz + bz 2 + uz 4 ^ 



where 



2\2 - "2°1 = 6w/&, 



(42) 



according to fl35l) . Three of scalar products can be calculated directly by means of the same 
technique as for the first Jordan cell. The result is: 

«*2,o|*2^)) = 7rc 2i0 iV 2i 2/A 3 ; 
((^2,il^2,i)) = ^2 2 i/A 3 ; 

«*2,o|*2,l» = °> 



in full agreement with ( fT71) : the normalization constants N 2> %, N 2l i can be chosen as necessary. 

Two other scalar products can be also brought in correspondence with (fTTj) . if the inte- 
gration constants are fixed suitably, and the relation ( 142|) is taken into account. The first of 
them, up to normalization constant: 



((#2,i 1*2,2)) ~ J exp (^-(Xzz + bz 2 +uz 4 )j 



1 1 2 2ag + ag 2a^ + b _ + 6) 

-(zzjz -z — 



2A 



2A 



A 2 



A 3 



dzdz\ 



(°) _j_ ^X 1 ) 



vanishes, if the constants satisfy the relation 1 + ct 2 
gives expressions for constants: 



-26, which together with ((42 



(o) 
"2,1 



3w (i) , 3w 
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The second scalar product is: 



((^2,2 1*2,2)) ~ j exp f-(Xzz + bz 2 + W2 4 )j 



X 1 ) , U\2 

4 - 1 A 4 v ~ 2,2 



6) 



(1)N2 



2:21 , _s2 1 a 2,2 + ^ 2 , 2/?2,2 / _%2 



a 



(i) 
2,2 



A 2 A 2 

2(ag + 6)ag _ 2w, _. 3 ' 

(zzjz — zz 7-{zz) 

A 



A 2 



dzdz. 



A A 3 

The condition, that ( H3|) vanishes, provides the expression for the integration constant: 

9u 2 



(43) 



4/3 : 



2,2 



fe 2 



+ 12w + b z 



and finally, the associated function fy 2 ,2 can be written in a rather compact form: 

2 



2.2 



Az + (6-^)z-2^ 3 ) + 18c(l-^ 2 + ^ 
b J b b 



exp 



■(Xzz + bz 2 + uz 4 ) . (44) 



Therefore, for the second Jordan cell all scalar products ( ITT)) are under control. Thus, we 
have elaborated the algorithm for construction, step by step, of higher Jordan cells with 
dimensionalities p n = n + 1 and energy E n = 2X(n + 1), although the explicit expressions for 
ty n k will become more and more complicated. 



Appendix. 

The set of relevant integrals over z, z was obtained in [16], and it'll be given again below for 
the reader convenience. The basic integral [30], [16] is: 

J(A, b, c) = J exp [-(Xzz + bz 2 + cz 2 )}dzdz = 2 7 r<r 1 ; 8 = ^(A 2 -46c), 

with A > (b + c) . Then, the next required integrals with power pre-exponential integrand 



arc: 



ijv.M = y z N z M exp [— (Azz + te 2 )]<iz<i2 



with X > b and integer AT, M. These integrals vanish for odd values of (N + M) due to 
antisymmetry under a space reflection (xi,x 2 ) — >■ — (xi,x 2 ). In turn, for even values of 
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(N + M) the integrals can be calculated by suitable differentiations of b, c) = 1(5) : 



with (L) P = L(L + 1)...(L + P - 1). In particular, it is clear that I N>M = for M > N. 
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